Measurements of Cross Sections and Charged Pion Spectra 
in proton— carbon Interactions at 31 GeV/c 
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Interaction cross sections and charged pion spectra in p+C interactions at 31 GeV/c were mea- 
sured with the large-acceptance NA61/SHINE spectrometer at the CERN SPS. These data are 
required to improve predictions of the neutrino flux for the T2K long-baseline neutrino oscillation 
experiment in Japan. A set of data collected during the first NA61/SHINE run in 2007 with an 
isotropic graphite target with a thickness of 4% of a nuclear interaction length was used for the 
analysis. The measured p+C inelastic and production cross sections are 257.2 ± 1.9 ± 8.9 and 
229.3 ± 1.9 ± 9.0 mb, respectively. Inclusive production cross sections for negatively and positively 
charged pious are presented as functions of laboratory momentum in ten intervals of the laboratory 
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polar angle covering the range from up to 420 mrad. The spectra are compared with predictions 
of several hadron production models. 

PACS numbers: 13.85.Lg,13.85.Hd,13.85.Ni 

Keywords: p+C interaction, inelastic cross section, inclusive pion spectra 



I. INTRODUCTION 

The NA61/SHINE (SPS heavy ion and neutrino ex- 
periment) experiment at the CERN Super Proton Syn- 
chrotron (SPS) pursues a rich physics program in vari- 
ous fields lUHI. First, precise hadron production mea- 
surements are performed for improving calculations of 
the neutrino flux in the T2K neutrino oscillation ex- 
periment [5j, as well as for more reliable simulations of 
cosmic-ray air showers in the Pierre Auger and KAS- 
CADE experiments [HllZ]- Second, p-l-p, p-l-Pb and nu- 
cleus-|-nucleus collisions will be studied extensively at 
SPS energies. A collision energy and system size scan 
will be performed with the aim to study properties of 
the onset of deconfinement and searching for the critical 
point of strongly interacting matter. This paper presents 
the first NA61/SHINE results on charged pion spectra 
in p+C interactions at 31 GeV/c, which are needed for 
an accurate neutrino flux prediction in the T2K experi- 
ment. The results are based on the data collected during 
the first run in 2007. NA61/SHINE is a large-acceptance 
hadron spectrometer [T]. Its main detector components, 
software, calibration, and analysis methods were inher- 
ited from the NA49 experiment [H |9j . 

T2K is a long-baseline neutrino experiment in Japan, 
which uses a high-intensity neutrino beam produced at 
J-PARC.^ It aims to precisely measure the t'^ — )■ 
appearance and disappearance [S'. In order to gen- 
erate the neutrino beam a high-intensity 30 GeV (ki- 
netic energy) proton beam impinging on a 90-cm-long 
graphite target is used, where tf and K mesons decaying 
into (anti)neutrinos are produced. The neutrino fluxes 
and spectra are then measured both at the near-detector 
complex, 280 m from the target, and by the Super- 
Kamiokande (SK) detector located 295 km away from 
the neutrino source and 2.5° off-axis. Neutrino oscilla- 
tions are probed by comparing the neutrino flux mea- 
sured at SK to the predicted one. In order to predict the 
flux at SK one uses the near-detector measurements and 
extrapolates them to SK with the help of Monte Carlo 
simulations. Up to now, these Monte Carlo predictions 
have been based on hadron production models only. For 
more precise predictions, which would allow the reduc- 
tion of systematic uncertainties to the level needed for 
the T2K physics goals, measurements of pion and kaon 
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Figure 1: (Color online) The prediction from the T2K beam 
simulation: the {p, 6} distribution for positively charged pi- 
ous weighted by the probability that their decay produces a 
muon neutrino passing through the SK detector. 



production off carbon targets are essential [TH3]. The 
purpose of the NA61/SIIINE measurements for T2K is to 
provide this information at exactly the proton extraction 
energy of the J-PARC Main Ring synchrotron, namely, 
30 GeV kinetic energy (approximately 31 GeV/ c momen- 
tum). Presently, the T2K neutrino beam- line is set up 
to focus positively charged hadrons, in such a way that 
it produces a beam. Spectra of positively charged 
pious presented in this paper constitute directly an es- 
sential ingredient in the neutrino flux calculation. The 
kinematic region of interest for positively charged pions 
whose daughter muon neutrinos pass through the SK de- 
tector is shown in Fig. [l] for the kinematic variables p - 
the momentum of a given particle and 9 - its polar angle 
in the laboratory frame. 

In addition to providing reference data for T2K, precise 
results on particle production in proton-carbon interac- 
tions furnish important input to improve hadronic gener- 
ators needed for the interpretation of air showers initiated 
by ultra-high-energy cosmic particles (see, e.g., [TDHI3])- 

The NA61/SIIINE data wfll also allow testing and im- 
provement of existing hadron production models in an 
intermediate energy region that is not well constrained 
by measurements at present. 

Experimental data on proton-nucleus interactions in 
the region of momentum of a few tens of GeV/c are rather 
scarce. Incident proton momenta from 3 up to 12 GeV/c 
were explored by the HARP Collaboration [T3] for a large 
variety of nuclear targets, including carbon |15j . Almost 
40 years ago an experiment was performed at CERN at 
24 GeV/c [TB] . Several targets were used but the momen- 
tum and angular range of produced particles was limited 
to angles from 17 to 127 mrad and momenta from 4 to 
18 GeV/ c. Recently the MIPP particle production exper- 
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iment at Fermilab presented preliminary results from its 
first data collected at 120 GeV/c [17]. Precise measure- 
ments of pion production in proton-carbon interactions 
at 158 GeV/c are available from the NA49 Collabora- 
tion 0. 

The paper is organized as follows. In Sec. II the 
NA61/SIIINE experimental set-up is described. Details 
on the beam, trigger, and event selection are given in 
Sec. III. Data reconstruction, simulation, and detector 
performance are described in Sec. IV. Analysis techniques 
and final results are presented in Sees. V and VI, respec- 
tively. These results are compared with hadron produc- 
tion models in Sec. VII. A summary in Sec. VIII closes 
the paper. 



II. THE NA61/SHINE SET-UP 

The NA61/SHINE experiment is a large-acceptance 
hadron spectrometer in the North Area H2 beam-line of 
the CERN SPS. The schematic layout is shown in Fig. [2] 
together with the overall dimensions. 

The main components of the current detector were con- 
structed and used by the NA49 Collaboration [B]. A set 
of scintillation and Cherenkov counters as well as beam 
position detectors (BPDs) upstream of the spectrome- 
ter provide timing reference, identification, and position 
measurements of the incoming beam particles. Details on 
this system are presented in Sec. |III| The main tracking 
devices of the spectrometer are large- volume time projec- 
tion chambers (TPCs). Two of them, the vertex TPCs 
(VTPC-1 and VTPC-2 in Fig.|2|, are located in a free gap 
of 100 cm between the upper and lower coils of the two 
superconducting dipole magnets. Their maximum com- 
bined bending power is 9 T-m. In order to optimize the 
acceptance of the detector at 31 GeV/c beam momen- 
tum, the magnetic field used during the 2007 data-taking 
period was set to a bending power of 1.14 T-m. Two 
large-volume main TPCs (MTPC-L and MTPC-R) are 
positioned downstream of the magnets symmetrically to 
the beam line. The TPCs are filled with Ar:C02 gas mix- 
tures in proportions 90:10 for VTPCs and 95:5 for MT- 
PCs. The particle identification capability of the TPCs 
based on measurements of the specific energy loss dE/dx 
is augmented by time-of-flight measurements {tof) using 
time-of-fiight (ToF) detectors. The ToF-L and ToF-R 
arrays of scintillator pixels have a time resolution of bet- 
ter than 90 ps [8]. Before the 2007 run the experiment 
was upgraded with a new forward time-of-fiight detec- 
tor (ToF-F) in order to extend the acceptance. The 
ToF-F consists of 64 scintillator bars with photomulti- 
plier (PMT) readout at both ends resulting in a time 
resolution of about 115 ps. The target under study is 
installed 80 cm in front of the VTPC-1. The results pre- 
sented here were obtained with an isotropic graphite tar- 
get of dimensions 2.5(W)x 2.5(11) x2(L) cm and with a 
density of p = 1.84 g/cm^. The target thickness along 
the beam is equivalent to about 4% of a nuclear interac- 



tion length (Ai). 



III. BEAM, TRIGGER AND DATA SAMPLES 

A 31 GeV/c secondary hadron beam is produced from 
400 GeV protons extracted from the SPS in slow extrac- 
tion mode. The beam is transported along the H2 beam- 
line toward the experiment. Collimators in the beam 
line are adjusted to get an average beam particle rate of 
15 kHz. The setup of beam detectors is illustrated in the 
inset on Fig. [2] Protons from the secondary hadron beam 
are identified by two Cherenkov counters, a CEDAR [18] 
and a threshold counter, labeled CI and C2, respectively. 
The CEDAR counter, using a six-fold coincidence, pro- 
vides positive identification of protons, while the thresh- 
old Cherenkov counter, operated at pressure lower than 
the proton threshold, is used in anti-coincidence in the 
trigger logic. The fraction of protons in the beam was 
about 14%. A selection based on signals from Cherenkov 
counters allowed the identification of beam protons with 
a purity of about 99%. A consistent value for the pu- 
rity was found by bending the beam into the TPCs with 
the full magnetic field and using the dE/dx identification 
method. 

Two scintillation counters SI and S2 provide beam 
definition, together with the two veto counters VO and 
VI with a 1-cm-diameter hole, which are coUimating the 
beam on the target. The SI counter provides also the 
timing (start time for all counters). Beam protons are 
then selected by the coincidence SI • S2 • V(5 • VT • CI • C2. 
The trajectory of individual beam particles is measured 
in a telescope of beam position detectors along the beam 
line (BPD-1, -2 and -3 in Fig. [2]). These counters are 
small (3x3 cm) proportional chambers with cathode strip 
readout, providing a resolution of about 200 fim in two 
orthogonal directions; see [5] for more details. The beam 
profile and divergence obtained from the BPD measure- 
ments are presented in Fig. |3] 

The beam momentum was measured directly in a ded- 
icated run by bending the incoming beam particles into 
the TPCs with the full magnetic field. The measured 
beam momentum distribution is shown in Fig. [4] The 
mean value of 30.75 GeV/c agrees with the set value of 
30.92 GeV/ c within the available precision of setting the 
beam magnet currents («0.5%) in the 112 beam-line. 

Interactions in the target are selected by an anti- 
coincidence of the incoming beam protons with a small, 
2-cm-diameter, scintillation counter (S4) placed on the 
beam trajectory between the two vertex magnets (see 
Fig. [2]) . This minimum bias trigger is based on the dis- 
appearance of the incident proton. A measurement of 
the interaction cross section is discussed in Sec. IV Bl 

The results presented in this paper are based on the 
analysis of 667 x 10^ proton interaction triggers recorded 
with the carbon target inserted and 46 x 10^ proton inter- 
action triggers recorded with the carbon target removed. 
Additionally 47 x 10"^ events were recorded with beam 
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Figure 2: (Color online) The layout of the NA61/SHINE experiment at the CERN SPS (top view, not to scale). The chosen 
right-handed coordinate system is shown on the plot. The incoming beam direction is along the z axis. The magnetic field 
bends charged particle trajectories in the x-z (horizontal) plane. The drift direction in the TPCs is along the y (vertical) axis. 



proton triggers and 33 x 10"^ with beam triggers without 
any beam particle identification. 



IV. DATA RECONSTRUCTION, SIMULATION 
AND DETECTOR PERFORMANCE 

A. Calibration 

The calibration procedure of the 2007 NA61 data was 
largely based on the approach developed for the NA49 
experiment [5] and consists of several consecutive steps 
resulting in optimized parameters for 

(i) detector geometry, drift velocity, and residual cor- 
rections, 

(ii) magnetic field, 

(iii) time-of-flight measurements, and 

(iv) specific energy loss measurements. 

Each step involved reconstruction of the data required 
to optimize a given set of calibration constants followed 
by verification procedures. Details of the procedure and 
quality assessment are presented in Ref. |4|. The qual- 
ity of detector calibration in quantities relevant for this 
paper is illustrated in the following subsections. 



B. Reconstruction 

Reconstruction algorithms used for the analysis de- 
scribed here are based on those developed by the NA49 
Collaboration [5] . The main steps of the reconstruction 
procedure are 

(i) cluster finding in the TPC raw data, and calculation 
of a cluster center of gravity and total charge, 

(ii) reconstruction of local track segments in each TPC 
separately, 

(iii) matching of track segments from different TPCs 
into global tracks, 

(iv) track fitting through the magnetic field and deter- 
mination of track parameters at the first measured 
TPC cluster, 

(v) determination of the interaction vertex as the in- 
tersection point of the incoming beam particle with 
the middle target plane, 

(vi) refitting the particle trajectory using the interaction 
vertex as an additional point and determining the 
particle momentum at the interaction vertex, and 

(vii) matching of ToF-F hits with the TPC tracks. 

An example of a reconstructed event is shown in Fig. [5] 
One can observe different track topologies, including long 
tracks hitting the ToF-F. 
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Figure 3: (Color online) Top: The beam spot as measured by 
BPD-3 after the VI cut described in the text. Bottom: The 
beam divergence in x and y. 



C. Monte Carlo simulation 

Simulation of the NA61 detector response used to cor- 
rect the raw data consists of the following steps (see [T5] 
for more details): 

(i) generation of p+C interactions at 31 GeV/ c using 
the VENUS4.12 Uni model, 

(ii) propagation of outgoing particles through the de- 
tector material using the GEANT 3.21 package [21], 
which takes into account the magnetic field as well 
as relevant physics processes, such as particle inter- 
actions and decays, 

(iii) simulation of the detector response using dedicated 
NA61 packages which introduce distortions corre- 
sponding to all corrections applied to the real data 
(see Sec.[!VAl), 



(iv) storage of the simulated events in a file which has 
the same format as the raw data, 

(v) reconstruction of the simulated events with the 
same reconstruction chain as used for the real data, 
and 

(vi) matching of the simulated and reconstructed tracks. 
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Figure 4: (Color online) The beam momentum distribu- 
tion measured by the reconstruction of beam particles in the 
TPCs. Only statistical errors are shown. 



Finally, the ratio of the number of reconstructed 
(matched) tracks in a given {p, 0) bin to the number 
of generated tracks of the specific particle type in this 
bin is used as the global correction factor (Sec. VG). Al- 
ternatively, flat-phase-space Monte Garlo simulation was 
used to calculate the same ratio for a given type of pri- 
mary particle and similar results were obtained within 
errors for the geometrical acceptance and reconstruc- 
tion efficiency corrections. The simulation based on the 
VENUS4.12 model, which has realistic ratios between the 
different particle species, was used to obtain the correc- 
tions due to the admixture of secondary particles e.g., 
from weak decays of other particles and from secondary 
interactions. The systematic uncertainties due to the de- 
pendence of the corrections on the model are discussed 
in Sec. IVGl 



D. Detector performance 



The quality of measurements was studied by recon- 



structing the momentum of beam particles (see Sec. Ill) 
and masses of and A particles from their V'^ decay 
topology. As an example the invariant mass distribution 
of K'g candidates is plotted in Fig. [6j The peak position 
(0.4976 ± 0.0005 GeV/c^) agrees with the known 
mass. The widths of the K'g and A invariant mass peaks 
are reasonably reproduced by the Monte Garlo simula- 
tion, which confirms its validity for the studies presented 
below. 

The track reconstruction efficiency and resolution of 
kinematic quantities were calculated by matching of sim- 
ulated and reconstructed tracks. As an example, the re- 
construction efficiency as a function of momentum for 
negatively charged tracks in the polar angle interval 
[100,140] mrad is shown in Fig. [7| The momentum reso- 
lution a{p)/p'^ , averaged over different track topologies in 
the detector, was estimated to be about 2x 10"'^, 7x10""^, 
and 3x10-2 (GeV/c)-i at p > 5GeV/c, p = 2GeV/c, 



Figure 5: (Color online) Topview (projection to the x-z plane) of a p+C interaction as detected by the four NA61 TPCs and 
the ToF-F detector (for details see Sec. and Fig. |2|. The red lines correspond to the trajectories of tracks produced in the 
interaction and reconstructed using the TPC clusters indicated by yellow points. Stars correspond to the reconstructed ToF-F 
hits. 



and p — 1 GeV/ c, respectively. The results depend some- 
what on particle production properties, and the quoted 
numbers refer to negatively charged particles (more than 
90% of them are pions) produced by Venus and passing 
the track selection cuts described in Sec. IV Al 

Particle identification was performed based on the en- 
ergy loss measurements in the TPCs and on the time- 
of-flight information from the ToF-F detector. The cali- 
brated dE/dx distributions as a function of particle mo- 
mentum for positively and negatively charged particles 
are presented in Fig. [8] The Bethe-Bloch parametriza- 
tion of the mean energy loss, scaled to the experimental 



data (see Sec. VD), is shown by the curves for positrons 
(electrons), pions, kaons, protons, and deuterons. The 
typical achieved dE/dx resolution is 3%-5%; see Fig. [9j 

The intrinsic tof resolution of the ToF-F detector is 
about 115 ps as derived from measurements of tof for 
particles traversing the overlap region of two adjacent 
scintillator bars (see Fig. 10). Figure [TT] shows distribu- 
tions of the mass squared rrr (derived from the to/, fitted 
momentum, and path length) versus dE/dx (measured 
in the TPCs) in selected intervals of particle momentum. 



The combined rn? and dE/dx measurements allow us to 
extract yields of identified particles even in the cross-over 
region of the Bethe-Bloch curves (1-4 GeV/c momentum 
range) . 

The geometrical acceptance of the ToF-F detector is 
limited to particle momenta above about 0.8 GeV/c. 
However, in the low-momentum region (less than about 
1 GeV/ c) the dE/dx information alone is sufficient to dis- 
tinguish pions from electrons/positrons, kaons, and pro- 
tons; see Fig. [8j 

About 47% of all charged particles produced at the in- 
teraction vertex pass the standard TPC track selection 
cuts (see Sec V A). In particular, the uninstrumented re- 



gion around the beam-line in VTPC-1 and VTPC-2 (see 
Fig. [2]) and the magnet apertures lead to a limited accep- 
tance in azimuthal angle. This is illustrated in Fig. [T2| 
where a distribution of all accepted negatively charged 
particles is shown as a function of azimuthal and polar 
angles. The fraction of charged particles accepted in the 
TPCs and in the TPCs plus ToF-F is plotted as a func- 
tion of momentum and polar angle in Fig. |13[ As the 
azimuthal acceptance varied with the polar angle 9 an 
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Figure 6: (Color online) Invariant mass distribution of recon- 
structed Kg candidates. The fitted mean value of the peak 
is given in the legend. The dashed histogram shows the sim- 
ulated distribution, which is normalized to the data on the 
right side of the peak. 
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Figure 7: (Color online) Track reconstruction efficiency for 
negatively charged particles as a function of momentum in 
the polar angle interval [100,140] mrad. 



individual wedge in azimuthal angle was adjusted for 
each particular 9 bin. Comparing plots in Fig. [13] with 
Fig. [l] one concludes that the NA61/SHINE acceptance 
fully covers the phase space region of interest for T2K. 



V. ANALYSIS TECHNIQUES 

This section presents the procedures used for data 
analysis. Crucial for this analysis is the identification 
of the pions produced. Depending on the momentum in- 
terval, different approaches were adopted, which led also 
to different track selection criteria. The task is facilitated 
for the negatively charged pions by the observation that 
more than 90% of primary negatively charged particles 
produced in p-l-C interactions at this energy are 7r~ , and 
thus the analysis of 7r~ spectra can also be carried out 
without additional particle identification. 

This section starts from a presentation of event and 




10 

p [GeV/c] 

Figure 8: (Color online) Specific energy loss in the TPCs for 
positively (top) and negatively {bottom) charged particles as 
a function of momentum. Curves show parametrizations of 
the mean dE/dx calculated for different particle species. 



track selection cuts. Next, the analysis of interaction 
cross sections is described [22] • Finally, three analysis 
methods applied to obtain pion spectra are introduced. 
These are 

(i) analysis of tt" mesons via measurements of neg- 
atively charged particles {h~ analysis |23|): see 
Sec.[VCl 

(ii) analysis of 7r+ and 7r~ mesons identified via dE/dx 
measurements in the TPCs {dE/dx analysis at low 
momentum [21]); see Sec. [VP] and 

(iii) analysis of tt"*" and 7r~ mesons identified via time-of- 
flight and dE/dx measurements in the ToF-F and 
TPCs, respectively {tof-dE/dx analysis [25j); see 
Sec.lVEl 

Each analysis yields fully corrected pion spectra with 
independently calculated statistical and systematic er- 
rors. The spectra were compared in overlapping phase- 
space domains to check their consistency. Complemen- 
tary domains were combined to reach maximum accep- 
tance. 

The final results refer to pions (denoted as primary 
pions) produced in p-|-C interactions at 31 GeV/c by 
strong interaction processes and in the electromagnetic 
decays of produced hadrons. 
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Figure 9: (Color online) The dE/dx distributions for posi- 
tively (full line) and negatively (dashed line) charged parti- 
cles at 5 and 15 GeV/ c particle momenta; the momentum bin 
width is 0.4 GeV/c. 



A. Event and track selection 

This section presents event and track selection criteria 
common for all analysis methods. Selection criteria spe- 
cific to a given method are described in the corresponding 
subsection belovif. 

The analysis was based on a sample of 521 x 10^ events 
selected from the total sample of 667 x 10^ registered 
and reconstructed proton interaction triggers recorded 
with the carbon target inserted. The selected events are 
required to have signals in each plane of all three BPD 
detectors with properly reconstructed beam tracks which 
include measured points on both planes of BPD-3. This 
criterion essentially removes contamination by interac- 
tions upstream of the target. For the event sample with 
the target removed, the selection reduces the number of 
events from 46 x 10^ to 17 x 10^. 

In order to select well-measured tracks in the TPCs as 
well as to reduce the contamination of tracks from sec- 
ondary interactions and weak decays the following track 
selection criteria were applied: 

(i) the track momentum fit at the interaction vertex 
should have converged, 

(ii) the total number of reconstructed points on a track 
should be at least 30, 
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Figure 10: (Golor online) Top: Distribution of the difference 
between a particle's time-of-flight measured independently by 
two overlapping scintillator bars of the ToF-F detector. The 
width of the distribution is about 160 ps, indicating a tof 
resolution of about 115 ps for a single measurement. Bottom: 
Mass squared, derived from the ToF-F measurement and the 
fitted path length and momentum, versus momentum p. The 
lines show the expected mass squared values for different par- 
ticles. 



(iii) the sum of the number of reconstructed points in 
VTPC-1 and VTPC-2 should be at least 12, 

(iv) the ratio of the total number of reconstructed points 
to the maximum possible number of points derived 
from the track trajectory with respect to the detec- 
tor geometry should be larger than 0.5, 

(v) the distance between a track extrapolated to the 
target plane and the interaction point (impact pa- 
rameter) should be smaller than 4 cm in both trans- 
verse directions separately, and 

(vi) the track azimuthal angle should be within an az- 
imuthal angle wedge around the horizontal plane; 
the wedge size depends on polar angle and the small- 
est one (±20° for /i^ and dE/dx analyses and ±10° 
for the tof-dE / dx analysis) was used for the largest 
polar angles; see Fig. [121 



B. Cross section measurements 

For normalization and cross section measurements we 
adopted the same procedure as the one developed by the 
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Figure 11: (Color online) Examples of two-dimensional m^-dE/dx plots for positively charged particles in three momentum 
intervals indicated in the panels. 2a contours around fitted pion peaks are shown. The left and middle plots correspond to the 
dE/dx cross-over region while the right plot is at such a high momentum that the ToF-F resolution becomes a limiting factor. 
The combination of both measurements provides close to 100% purity in the pion selection over the whole momentum range. 
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Figure 12: (Color online) Polar angle (6) vs azimuthal angle 
(0) distribution for reconstructed negatively charged particles 
in the momentum interval 0.5 < p [GeV/c] < 5.0. 



NA49 Collaboration The minimum bias trigger on 
proton interactions, described in Sec. 



Ill allows us to 



define a "trigger" cross section which is used both for 
the normalization of the differential inclusive pion dis- 
tributions and for the determination of the inelastic and 
production cross sections. 

From the numbers of selected interactions, fulfilling the 
requirements on BPD signals and reconstruction of the 
proton beam particles as detailed in Sec. |V A[ we com- 
pute an interaction probability of (6.022 ± 0.034)% with 
the carbon target inserted and of (0.709 ± 0.007)% with 
the carbon target removed. These measurements lead 
to an interaction probability of (5.351 ± 0.035)% in the 
carbon target, taking into account the reduction of the 
beam intensity in the material along its trajectory. The 
corresponding "trigger" cross section is (298.1 ± 1.9 ± 
7.3) mb, after correcting for the exponential beam at- 
tenuation in the target. The systematic error on the 
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Figure 13: (Color online) Fraction of accepted particles as a 
function of momentum and polar angle, after the track accep- 
tance cuts (see Sec. VAl for negatively charged tracks (top), 
and after an additional ToF-F acceptance cut (see Sec. |VE[ ) 
for positively charged tracks (bottom). The first polar angle 
bin, [0,20] mrad, is fully covered by accepted particles up to 
7.6 GeV/c. 



"trigger" cross section was conservatively evaluated by 
comparing this value with the one obtained without any 
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event selection criteria. 

Two classes of processes are considered, inelastic and 
production interactions. The production processes are 
defined as those in which new hadrons are produced. The 
inelastic processes include in addition interactions which 
result only in disintegration of the target nucleus (quasi- 
elastic interactions). 

The inelastic cross section ainei was defined as the sum 
of all processes due to strong p+C interactions except 
coherent nuclear elastic scattering. Thus it includes in- 
teractions with production of new hadrons (production 
processes) and quasi-elastic interactions which lead only 
to break up of the carbon nucleus. The inelastic cross 
section was derived from the "trigger" cross section by 
applying two corrections: 

(i) subtraction of the contribution from coherent elastic 
scattering (47.2 ± 0.2 ± 5.0 mb), i.e., removal of 
those events in which the incoming beam particle 
undergoes a large-angle coherent elastic scatter on 
the carbon nuclei and misses the S4 counter, and 

(ii) addition of the lost inelastic events because 
of emitted charged particles hitting the S4 
counter (5.7 ± 0.2 ± 0.5 mb for protons and 
0.57 ± 0.02 ± 0.35 mb for pions and kaons). 

The corrections were calculated based on the 
GEANT4 [27^ simulation (with the QGSP_BERT physics 
list) of the beam line setup using the measured pro- 
file and divergence of the proton beam. The GEANT4 
angular distributions for coherent elastic scattering and 
quasi-elastic scattering, as well as the total inelastic cross 
section, were cross-checked against available experimen- 
tal measurements. The largest discrepancy of 7.5% was 
found for the coherent elastic events in which the scat- 
tered proton is outside of the S4 acceptance. This dis- 
crepancy was taken into account in the calculation of 
the systematic error, which includes also uncertainties in 
other simulation parameters, namely, beam position and 
divergence, as well as the size of the S4 counter. The total 
systematic error on the inelastic cross section amounts to 
8.9 mb. 

The production cross section is determined from the 
inelastic cross section by subtracting the cross section 
of quasi-elastic p-l-C interactions at 31 GeV/c which 
amounts to 27.9 ± 1.5 (sys) mb according to Glauber 
model calculations [28 . 

Details of the cross section analysis procedure can be 
found in Ref. [221 ■ 



C. The h analysis 

More than 90% of primary negatively charged parti- 
cles produced in p-l-C interactions at 31 GeV/c are neg- 
atively charged pions. Thus 7r~ meson spectra can be 
obtained by subtracting the estimated non-pion contribu- 
tion from the spectra of negatively charged particles and 



additional particle identification is not required. Note 
that this method is not applicable to the analysis of tt"*" 
meson spectra because of the much larger and unknown 
contribution of protons and K'^ mesons to all positively 
charged particles. 
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Figure 14: (Color online) The relative contribution of ac- 
cepted primary K~ and p to accepted primary negatively 
charged pions in different (p, 9) intervals calculated within 
the Venus model. 

First, spectra of all accepted negatively charged par- 
ticles were obtained in (p, 6) bins. The event and track 



selection criteria presented in Sec. VA were applied. 

Second, the Monte Carlo simulation described in 
Sec. IIV CI was used to calculate corrections for the con- 
tribution of electrons and primary K~ and p as well as 
secondary particles from weak decays (feed down), sec- 
ondary interactions, and photon conversions in the target 
and the detector material. Furthermore, the corrections 
include track reconstruction efficiency and resolution as 
well as losses due to the limited geometrical acceptance 
of the detector. Bin-by-bin correction factors were calcu- 
lated as the ratio of all generated primary tt~ mesons to 
all reconstructed and accepted negatively charged parti- 
cles in a given bin. 

Different contributions to the correction factors were 
studied separately. The losses due to the track recon- 
struction efficiency and the limited geometrical detector 
acceptance are shown in Figs. [T] and 13 respectively. The 
relative contribution of accepted primary and p to 
accepted primary negatively charged pions is presented 
in Fig. [T4j Finally, the relative contribution of accepted 
electrons and non-primary negatively charged hadrons 
to accepted primary negatively charged pions is given in 
Fig. [15] The correction for secondary pions and electrons 
is largest for small momenta. The former are mainly due 
to strange particle decays close to the primary vertex; 
the latter originate mostly from conversion of tt" decay 
photons. Uncertainty due to this correction is discussed 
in Sec. IVGl 

The inverse correction factor and the geometrical ac- 
ceptance versus momentum for the polar angle interval 
[140,180] mrad are shown in Fig. 16 as an example. 
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Figure 15: (Color online) The relative contribution of ac- 
cepted electrons and non-primary negatively charged hadrons 
to accepted primary negatively charged pions in different 
(p, 6) intervals calculated within the Venus model. The low- 
momentum region where the ratio is larger than 3 is excluded. 
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Figure 16: (Color online) Example of momentum depen- 
dence of the inverse correction factor and geometrical accep- 
tance (top) and systematic errors {bottom) for the anal- 
ysis for negatively charged pions in the polar angle interval 
[140,180] mrad. For details concerning systematic uncertain- 
ties, see Sec. |V G[ 



D. The dE/dx analysis at low momentum 

The analysis of charged pion production at low mo- 
mentum was done by means of particle identification via 
measurements of specific energy loss in the TPCs. Mea- 
surements of tof are not available for low-momentum 
particles since they do not reach the ToF-F detectors. 
A reliable identification of tt"*" mesons was not possible 
at momenta above 1 GeV/ c where the Bethe-Bloch (BB) 
curves for pions, kaons, and protons cross each other (see 
Fig. [8|. On the other hand, for tt" mesons, where the 
contribution of and antiprotons is almost negligible, 
the dE/dx analysis could be extended in momentum up 
to 3 GeV/c allowing consistency checks with the other 
analysis methods in the region of overlap. 

The procedure of particle identification, described be- 
low, is tailored to the region where a fast change of en- 
ergy loss with momentum is observed. In order to op- 
timize the parametrization of the BB function, samples 
of e*, TT^, ^ p, and d tracks with reliable particle 
identification were chosen in the /3"f range from 0.2 up 
to 100. The dependence of the BB function on f3'y was 
then fitted to the data using the Sternheimer and Peierls 
parametrization of Ref. |29j . This function was then used 
to calculate for every track of a given momentum the 
expected [dE /dx)BB values for all possible identifica- 
tion hypotheses for comparison with the measured mean 
{dE / dx)iiata- A small (a few percent) dependence of the 
mean {dE/dx) data on the track polar angle had to be 
corrected for. 

The identification procedure was performed in (p, 9) 
bins. Narrow momentum intervals (of 0.1 GeV/c for 
p < 1 GeV/c and 0.2 GeV/c for 1 < p < 3 GeV/c) were 
chosen to account for the strong dependence of dE/dx on 
momentum. The event and track selection criteria pre- 



sented in Sec. V A were applied. In each {p, 9) bin an un- 
binned maximum likelihood fit (for details see Ref. [50] ) 
was performed to extract yields of tt"*" and tt" mesons. 
The probability density functions were assumed to be a 
sum of Gaussian functions for each particle species, cen- 
tered on {dE/dx)BB with variances derived from data. 
The dE/dx resolution is a function of the number of mea- 
sured points and the particle momentum. In the 7r+ anal- 
ysis three independent abundances were fitted (7r+, 
and proton) while in the tt" analysis we were left with 
only two independent abundances (tt" and K~). The 
e"*" and e~ abundances were determined from the total 
number of particles in the fit. 

As an example, the dE/dx distribution for positively 
charged particles in the momentum bin [0.7,0.8] GeV/c 
and angular bin [180,240] mrad is shown in Fig. 17 com- 
pared with the distribution obtained using the fitted pa- 
rameters. 

Finally, the Monte Carlo simulation described in 
Sec. |IVC| was used to calculate bin-by-bin corrections 
for pions from weak decays and interactions in the tar- 
get and the detector material. The corrections include 
also track reconstruction efficiency and resolution as well 
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as losses due to pion decays and the limited geometrical 
acceptance of the detector. The inverse correction factor 
and geometrical acceptance versus momentum for the po 
lar angle interval [140,180] mrad is shown in Fig 
an example. 



individual as well as combined global correction factors 
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as 



E. The tof-dE/dx analysis 

High purity particle identification can be performed by 
combining the tof and dE/dx information. Moreover, in 
the momentum range 1-4 GeV/c, where dE/dx bands 
for different particle species overlap, particle identifica- 
tion is in general only possible using the tof method; see 
Figs. [8] and 11 for illustration. The ToF-F detector was 
designed to cover the necessary acceptance in momentum 
and polar angle required by the T2K experiment. 

First, a particle mass squared, m? , was calculated us- 
ing the tof , momentum, and path length measurements. 
The m^-dE/dx distributions of all accepted positively 
(and negatively) charged particles were then obtained in 
each (p, 9) bin. The event and track selection criteria 
presented in Sec. V A were applied. In addition only 
tracks which gave a signal in the ToF-F detector were 
selected. Two independent topologies of tracks, emitted 
to the left and to the right with respect to the incoming 
beam direction, have different acceptances. The analysis 
was performed for each topology separately. When both 
topologies are present the final result is the weighted av- 
erage of the two. Momentum-dependent azimuthal angle 
cuts were applied such that in these angular intervals the 
detector acceptance was close to 100%. These cuts select 
tracks with a large number of measured points as well, 
thus with a high reconstruction efficiency. The averaged 
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Figure 17: (Color online) The dE/dx distribution for posi- 
tively charged particles in the momentum bin [0.7,0.8] GeV/c 
and angular bin [180,240] mrad compared with the distribu- 
tion calculated using the fitted relative abundances. A small 
systematic shift between the fit curves (Gaussian peaks cen- 
tered on the Bethe-Bloch values) and data is observed. This 
effect results in a major contribution to the total systematic 
error in the particle identification procedure. 



are shown in Fig. 19 for illustration. 

The pion accumulations in these distributions were 
parametrized by a product of Gaussian functions in 
and dE/dx. In each (p, 6) bin the bin- by-bin maximum 
likelihood method was applied to fit yields of 7r+ and tt^ 
mesons. In the fit the signal shapes parametrized in the 
first step were allowed to vary within narrow limits. The 
pion yields were calculated summing all particles within 
2(7 around the fitted pion peak; see Fig. 11 In addition 
the integral of the fitted function was computed and used 
for the evaluation of the systematic uncertainty. 

Finally, the Monte Carlo simulation described in 
Sec. |IV C] was used to calculate corrections for pions from 
weak decays and interactions in the detector material and 
target, track reconstruction efficiency and resolution, and 
losses due to pion decays. The ToF-F detection efficiency 
was estimated by requiring that a track traversing the 
ToF-F wall generates a hit in the ToF-F. The ToF-F in- 
efficiency is due also to double hits and readout inefficien- 
cies. The bin-by-bin correction factors for different biases 
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Figure 18: (Color online) Example of the momentum depen- 
dence of the inverse correction factor and geometrical accep- 
tance (top) and systematic errors (bottom) for the dE/dx anal- 
ysis for positively charged pions in the polar angle interval 
[140,180] mrad. The small differences between geometrical 
acceptance and global correction are due to the fact that pion 
loss and feed down corrections, which are the most important 
ingredients of the global correction, cancel to a large extent. 
For details concerning systematic uncertainties see Sec. |V G[ 
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as an 



were calculated separately and applied to the data. The 
inverse correction factors versus momentum, total cor- 
rection factor and geometrical acceptance for the polar 
angle interval [40,60] mrad are shown in Fig. 
example. 



F. Derivation of spectra 



than the data statistics and the total statistical errors 



The procedures presented in Sees. VC VD[ and |V E| 



were used to analyze events with the carbon target in- 
serted (I) and with the carbon target removed (R). The 
corresponding corrected numbers of tt~ and 7r+ mesons 
in p bins and 9 intervals are denoted as An^ and An^, 
where a stands for tt~ and 7r+. Note that the same 
event and track selection criteria as well as corrections 
discussed in Sees. |V CI IVDl and IVEI were used in the 



analysis of events with the target inserted and removed. 
The latter events allow us to correct the results for the 
contribution of out-of-target interactions. 

Then, the differential inclusive cross section of tt"*" and 
7r~ mesons is calculated as 



daa 
dp 



f^tr 



1 Ani 



Ani 



1-e V^^-^ Ap iV« Ap 



(1) 



where 

(i) atrig = (298.1 ± 1.9 ± 7.3) mb is the "trigger" cross 
section as given in Sec. |VB[ 

(ii) and are the numbers of events selected (see 
Sec. VA) for the analysis of events with the target 



inserted and removed, respectively, 

(iii) Ap is the bin size in momentum, and 

(iv) e = 0.118±0.001 is the ratio of the interaction prob- 
abilities for operation with the target removed and 
inserted. 

The correction for the contribution of particles from 
out-of-target events (the term e/N^ ■ An^/ Ap in Eq. [l]) 
amounts on average to about 7% and 3% in the first two 
polar angle intervals. It is smaller than about 2% for 
polar angle bins above 40 mrad. 

The pion spectra normalized to the mean pion multi- 
plicity in production interactions was calculated as 



dria 
dp 



"^prod 



d(Ja 

dp 



(2) 



where aprod is the cross section for production processes. 



G. Statistical and systematic errors 

Statistical errors on the pion spectra include contribu- 
tions from the finite statistics of data and from the Monte 
Carlo simulation used to obtain the correction factors. 
The Monte Carlo statistics was about ten times larger 
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Figure 19: (Color online) Example of momentum dependence 
of the inverse correction factors (top), total inverse correc- 
tion factor and geometrical acceptance {middle), and system- 
atic errors {bottom) for the tof~dE/dx analysis for positively 
charged pious in the polar angle interval [40,60] mrad. e^ec 
and eToF are the efficiencies of the reconstruction and of the 
ToF-F, respectively. The feed-down correction accounts for 
pions from weak decays which are reconstructed as primary 
particles, while the pion loss accounts for pions lost due to de- 
cays or secondary interactions. For details concerning system- 
atic uncertainties, see Sec. VG For the tof-dE / dx analysis, 



within a given theta bin, in addition to the cut on the az- 
imuthal angle, a momentum-dependent selection on the track 
topology was applied. The purpose of this cut is to select 
only regions with a flat acceptance along the whole momen- 
tum range. 



13 



are dominated by the statistical uncertainty of the data. 

Systematic errors on the pion spectra were estimated 
by varying track selection and identification criteria as 
well as parameters used to calculate the corrections. The 
following track selection criteria were varied: the mini- 
mum number of points measured on the track, the az- 
imuthal angle, and the impact parameter cuts. It was 
found that the influence of such changes is small as com- 
pared to the statistical errors. The accuracy of correc- 
tions for acceptance and reconstruction efficiency was 
checked by comparison of the results obtained with inde- 
pendent track topologies as well as using two different al- 
gorithms for merging track segments from different TPCs 
into global tracks (track merging algorithm). The dom- 
inant contributions to the systematic error for all three 
analysis methods come from the uncertainty in the cor- 
rection for secondary interactions and for weak decays of 
strange particles. We assigned an uncertainty of 30% of 
the correction for both sources. The systematic error due 
to the admixture of pions from the decays of strange par- 
ticles reconstructed at the primary vertex depends on the 
knowledge of strange particle production. The estimate 
of this error was based on the following. 

(i) Comparison of the number of V'^ decays recon- 
structed in the data and in the Venus4.12 model 
with default parameters. 

(ii) Variation of the strange particle yields in dif- 
ferent Monte Carlo generators; for example the 
K~/tt^ ratio in p-l-C interactions at 31 GeV/c 
from FLUKA2008 [31], Urqmd1.3.1 [32] and 
GiBUU 1.3.0 [531 (with default value of physics pa- 
rameters) is 3.73%, 3.67% and 3.35%, respectively; 
these values can be compared with 4.06% from the 
Venus4.12 [5D] generator used in this paper for cal- 
culation of corrections for pions from weak decays. 

(iii) Comparison of the Venus4.12 predictions with the 
measured /tt' and K~^/n^ ratios from P-I-B4C 
interactions at 24 GeV/c at large momenta and 
small angles. 

In the h" analysis a 20% uncertainty in the correction 
for electron admixture to the pion sample was assumed. 
Examples of systematic errors for all three analyses are 
presented in the bottom panels of Figs. [16] [18] and [19] 
Losses of inelastic interactions because an emitted par- 
ticle hits the S4 counter were estimated to be negligibly 
small, namely, about 1% or smaller in the relevant phase 
space for the analysis. The total systematic error was 
calculated as a sum of different contributions added in 
quadrature. It does not include the overall uncertainty 
due to the normalization procedure, namely, 2.5% and 
2.3% for the normalization to the inclusive cross sec- 
tion and mean pion multiplicity in production events, re- 
spectively. Statistical and systematic uncertainties were 
added in quadrature in order to calculate the total error. 
The first quoted error always refers to the statistical and 
the second to the systematic uncertainty. 



VI. RESULTS 

This section presents results on inelastic and produc- 
tion cross sections as well as on differential spectra of 7r+ 
and TT~ mesons in p-l-C interactions at 31 GeV/c. 



A. Inelastic and production cross sections 

The corrections to the "trigger" cross section described 
in Sec. IV Bl result in a total inelastic cross section of 

cr,„e, = 257.2 ± 1.9 ±8.9 mb. 

The production cross section was calculated from the 
inelastic cross section by subtracting the quasi-elastic 
contribution. The result is 



^prod 



229.3 ± 1.9 ±9.0 mb. 



The production cross section is compared to previous 

The NA61/SHINE result at 



measurements in Fig. 20 



31 GeV/c is consistent with the measurement of Ref. [M] 
at 20 GeV/ c, after subtraction of the quasi-elastic contri- 
bution of 30.4 ±1.9 (sys) mb, and with that of Ref. [35] 
at 60 GeV/ c. The measurements presented in Ref. [55] at 
different momenta are found to be systematically larger, 
as was already noted in Ref. |35j , where the authors com- 
mented "We know of no reason for this discrepancy" . 
Note that the data from Ref. [5S] are about 30 mb higher 
than other results; if interpreted as measurements of in- 
elastic cross section, they would agree with the other pro- 
duction cross section data, after the correction for the 
quasi-elastic contribution. 



B. Spectra of 7r+ and tt mesons 

The 7r+ and tt~ spectra presented in this section re- 
fer to pions produced in strong and electromagnetic pro- 
cesses in p+C interactions at 31 GeV/c. 

The spectra normalized to the inclusive cross section 
are shown in Figs . [21] and [22] for negatively and positively 
charged pions, respectively. The spectra are presented as 
functions of particle momentum in ten intervals of the 
polar angle. Both quantities are calculated in the lab- 
oratory system. The chosen binning takes into account 
the available statistics of the 2007 data sample, detector 
acceptance, and particle production kinematics. 

The spectra obtained by different methods (described 
in Sec. |v]) are presented separately. The agreement be- 
tween them is, in general, better than 10%. Note that 
data points in the same (p, 9) bin from different analysis 
methods are statistically correlated as they result from 
the analysis of the same data set. In order to obtain 
the final spectra consisting of statistically uncorrelated 
points, the measurement with the smallest total error 
was selected. The corresponding numerical values are 
presented in the Appendix, Table [l] and are also available 
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The results are presented in Figs. [23] and [24] for the 
spectra of 7t~ and tt"*" respectively. As already presented 
in [30] , Gheisha simulations qualitatively fail to describe 
the NA61/SHINE measurements at all production angles 
and momenta (see also, e.g., [IH])- The Urqmd1.3.1 
model qualitatively disagrees with the data only at low 
momenta (p < 3 GeV/c) and polar angles below about 
140 mrad. The Venus4.12 and Fluka2008 models fol- 
low the data trend in all measured polar angle intervals. 



VIII. SUMMARY 



Figure 20: (Color online) Beam momentum dependence of 
the production cross section for p+C interactions. The 
NA61/SHINE result (filled square) is compared with previ- 
ous measurements: Bellettini et al. (circle) [34) . Carroll et 
al. (triangle) 05] and Denisov et al. (open squares) [36j . 
For the NAG 1/ SHINE point, the error bar indicates statisti- 
cal and systematic uncertainties added in quadrature. The 
result from Ref. [34] was recalculated by subtracting from the 
measured inelastic cross section a quasi-elastic contribution 
at 20 GeV/c of 30.4 ± 1.9 (sys) mb. 



from |37| , where the different contributions to systematic 
uncertainties (see Sec. V G ) are given s epa rate ly. 

The final spectra are plotted in Figs. 23 and 24 for tt^ 
and TT"*", respectively, while the corresponding fractional 
errors are presented in Figs. [25] and [26] For the purpose 
of a comparison of the data with model predictions, the 
spectra were normalized to the mean tt* multiplicity in 
all production interactions. This avoids uncertainties due 
to the different treatment of quasi-elastic interactions in 
models as well as problems due to the absence of predic- 
tions for inclusive cross sections. 

The ratio of the final spectra of 7r+ and tt~ is presented 

The 77"*" to 7r~ 



as a function of momentum in Fig. 27 



ratio is close to 1 at low momenta and increases with 
increasing momentum. 

VII. COMPARISON TO MODEL PREDICTIONS 

As a first application of the measurement presented in 
this paper, it is interesting to compare the 7r~ and ir'^ 
spectra in p-l-C interactions at 31 GeV/c to the predic- 
tions of event generators of hadronic interactions. Models 
that have been frequently used for the interpretation of 
cosmic ray data, i.e., Venus4.12 j20j, Fluka2008 (31j . 
Urqmd 1.3.1 [32 , and Gheisha2002 [31] were selected. 
They are part of the CORSIKA [39| framework for the 
simulation of air showers and are typically used to gen- 
erate hadron-air interactions at energies below 80 GeV. 
In order to assure that all relevant settings of the gen- 
erators are identical to the ones used in air shower sim- 
ulations, p-|-C interactions at 31 GeV/c were simulated 
within CORSIKA in the so-called interaction test mode. 



This work presents inelastic and production cross sec- 
tions as well as positively and negatively charged pion 
spectra in p-l-C interactions at 31 GeV/c. These data 
are essential for precise predictions of the neutrino flux 
for the T2K long-baseline neutrino oscillation experi- 
ment in Japan. Furthermore, they provide important 
input to improve hadron production models needed for 
the interpretation of air showers initiated by ultra-high- 
encrgy cosmic particles. The measurements were per- 
formed with the large-acceptance NA61/SHINE spec- 
trometer at the CERN SPS. A set of data collected 
with a 4% Ai isotropic graphite target during the pilot 
NA61/SHINE run in 2007 was used for the analysis. The 
p-l-C inelastic and production cross sections were found 
to be 257.2 ± 1.9 ± 8.9 and 229.3 ± 1.9 ± 9.0 mb, re- 
spectively. Negatively and positively charged pion spec- 
tra as a function of laboratory momentum in ten inter- 
vals of the polar angle were obtained using three differ- 
ent analysis techniques. The final spectra were compared 
with predictions of hadron production models. 

The data presented in this paper already provide im- 
portant information for calculating the T2K neutrino 
flux. Meanwhile, a much larger data set with both the 
thin (4% Ai) and the T2K rephca carbon targets was 
recorded in 2009 and 2010 and is presently being ana- 
lyzed. This will lead to results of higher precision for pi- 
ous and extend the measurements to other hadron species 
such as charged kaons, protons, Kg and A. The new 
data will allow a further significant reduction of the un- 
certainties in the prediction of the neutrino flux in the 
T2K experiment. 
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are presented as a function of laboratory momentum (p) in different intervals of polar angle (d). Results obtained using three 
analysis methods are presented by different symbols: blue open circles, h" analysis; red open squares, dE/dx analysis; and 
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Figure 22: (Color online) DifTerential cross sections for tt^ meson production in p+C interactions at 31 GeV/c. The spectra 
are presented as a function of laboratory momentum (p) in different intervals of polar angle (6). Results obtained using two 
analysis methods are presented by different symbols: red open squares, dE/dx analysis, and black full triangles, tof-dE/dx 
analysis. Error bars indicate only statistical uncertainties. 
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Figure 23: (Color online) Laboratory momentum distributions of tt~ mesons produced in production p+C interactions at 
31 GeV/c in different intervals of polar angle (9). The spectra are normalized to the mean 7r~ multiplicity in all produc- 
tion p+C interactions. Error bars indicate statistical and systematic uncertainties added in quadrature. The overall uncer- 
tainty (2.3%) due to the normalization procedure is not shown. Predictions of hadron production models, Fluka2008 (solid 
line), UrqmdI.3.1 (dashed line), and Venus4.12 (dotted line), are also indicated. 
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Figure 24: (Color online) Laboratory momentum distributions of vr"*" mesons produced in production p+C interactions at 
31 GeV/c in different intervals of polar angle (9). The spectra are normalized to the mean n'^ multiplicity in all produc- 
tion p+C interactions. Error bars indicate statistical and systematic uncertainties added in quadrature. The overall uncer- 
tainty (2.3%) due to the normalization procedure is not shown. Predictions of hadron production models, Fluka2008 (solid 
line), UrqmdI.3.1 (dashed line), and Venus4.12 (dotted line) are also indicated. 
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Figure 25; (Color online) Relative uncertainties on the n~ cross sections shown in the same binning as Fig. 23 Statistical 
(dotted line), systematic (dashed line), and total (solid line) uncertainties are indicated. The overall uncertainty due to the 
normalization procedure is not shown. 
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Figure 26; (Color online) Relative uncertainties on the tt cross sections shown in the same binning as Fig. 24 Statistical 
(dotted line), systematic (dashed line), and total (solid line) uncertainties are indicated. The overall uncertainty due to the 
normalization procedure is not shown. 
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Figure 27: (Color online) The momentum dependence of the vr"'" to tt~ ratio. Errors are calculated taking into account only 
statistical uncertainties of the spectra plotted in Figs. |23| and |24[ Predictions of hadron production models, Fluka2008 (solid 
line), Urqmd1.3.1 (dashed line), and Venus4.12 (dotted linejare also indicated. 
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Appendix 



Table I: The NA61/SHINE results for the differential 7r+ and tt" pro- 
duction cross section in the laboratory system, da'^ /dp, for p+C inter- 
actions at 31 GoV/c. Each row refers to a different {^piow < P < Pup, 
Oiow < < 9up) bin, where p and 8 are the pion momentum and polar 
angle in the laboratory frame, respectively. The central value as well 
as the statistical (A^tat) and systematic (Asys) errors are given. The 
overall uncertainty (2.5%) due to the normalization procedure is not 
included. 
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2.8 3.0 5.80 0.77 0.36 2.8 3.0 5.74 0.20 0.18 

3.0 3.2 6.04 0.78 0.38 3.0 3.2 5.24 0.19 0.24 

3.2 3.4 3.76 0.63 0.24 3.2 3.4 4.48 0.18 0.20 

3.4 3.6 4.31 0.66 0.28 3.4 3.6 3.80 0.16 0.17 

3.6 3.8 3.58 0.62 0.23 3.6 3.8 3.24 0.15 0.15 

3.8 4.0 3.05 0.57 0.19 3.8 4.0 3.02 0.15 0.14 

4.0 4.2 3.32 0.60 0.21 4.0 4.2 2.40 0.13 0.10 

4.2 4.4 3.37 0.62 0.22 4.2 4.4 2.20 0.13 0.10 

4.4 4.6 1.28 0.41 0.08 4.4 4.6 1.92 0.12 0.08 

4.6 4.8 2.05 0.49 0.13 4.6 4.8 1.58 0.11 0.07 

4.8 5.0 1.34 0.41 0.09 4.8 5.0 1.50 0.10 0.06 

5.0 5.2 1.20 0.39 0.08 5.0 5.2 1.141 0.093 0.050 

5.2 5.6 1.16 0.18 0.08 5.2 5.6 0.748 0.051 0.033 

5.6 6.0 0.83 0.16 0.05 5.6 6.0 0.661 0.049 0.028 

6.0 6.4 0.407 0.037 0.017 

6.4 6.8 0.314 0.035 0.013 

6.8 7.2 0.283 0.031 0.012 

7.2 7.6 0.223 0.030 0.009 

180 240 0.2 0.3 6.73 0.38 0.48 0.2 0.3 5.61 0.31 0.71 

0.3 0.4 10.66 0.43 0.58 0.3 0.4 9.3 0.4 1.0 

0.4 0.5 15.58 0.50 0.75 0.4 0.5 13.1 0.4 1.2 

0.5 0.6 17.33 0.52 0.77 0.5 0.6 15.7 0.5 1.2 

0.6 0.7 20.03 0.54 0.84 0.6 0.7 17.7 0.5 1.2 

0.7 0.8 21.01 0.61 0.89 0.7 0.8 18.8 0.5 1.1 

0.8 0.9 22.4 0.7 1.0 0.8 0.9 19.0 0.5 1.1 

0.9 1.0 22.1 0.6 1.0 0.9 1.0 20.0 0.5 1.1 

1.0 1.2 23.0 1.6 1.7 1.0 1.2 19.3 0.4 1.1 

1.2 1.4 21.8 1.5 1.6 1.2 1.4 18.17 0.35 0.97 

1.4 1.6 18.5 1.3 1.3 1.4 1.6 16.74 0.34 0.74 

1.6 1.8 15.8 1.3 1.1 1.6 1.8 14.25 0.32 0.60 

1.8 2.0 15.2 1.2 1.0 1.8 2.0 12.50 0.29 0.50 

2.0 2.2 12.5 1.1 0.8 2.0 2.2 10.46 0.27 0.42 

2.2 2.4 11.6 1.1 0.8 2.2 2.4 9.29 0.26 0.35 

2.4 2.6 7.48 0.86 0.49 2.4 2.6 7.80 0.23 0.29 

2.6 2.8 7.62 0.86 0.49 2.6 2.8 6.63 0.21 0.24 

2.8 3.0 6.51 0.81 0.44 2.8 3.0 5.63 0.20 0.21 

3.0 3.2 5.89 0.75 0.39 3.0 3.2 4.30 0.20 0.22 

3.2 3.4 3.83 0.63 0.24 3.2 3.4 3.47 0.18 0.17 

3.4 3.6 4.95 0.66 0.33 3.4 3.6 2.80 0.16 0.14 

3.6 3.8 3.54 0.62 0.25 3.6 3.8 2.36 0.15 0.12 

3.8 4.0 2.91 0.56 0.19 3.8 4.0 1.93 0.14 0.10 

4.0 4.2 2.31 0.52 0.15 4.0 4.2 1.49 0.12 0.08 

4.2 4.4 1.74 0.46 0.12 4.2 4.4 1.49 0.12 0.07 

4.4 4.6 1.28 0.41 0.09 4.4 4.6 1.26 0.11 0.07 

4.6 4.8 1.72 0.45 0.11 4.6 4.8 1.04 0.10 0.05 

4.8 5.0 0.826 0.092 0.044 

5.0 5.2 0.593 0.076 0.032 

5.2 5.6 0.471 0.047 0.025 

5.6 6.0 0.275 0.037 0.017 

6.0 6.4 0.181 0.033 0.010 

240 300 0.2 0.3 8.49 0.44 0.54 0.2 0.3 7.36 0.37 0.87 

0.3 0.4 13.98 0.53 0.66 0.3 0.4 12.4 0.4 1.2 

0.4 0.5 17.89 0.59 0.78 0.4 0.5 15.8 0.5 1.4 

0.5 0.6 20.21 0.61 0.84 0.5 0.6 17.7 0.5 1.3 

0.6 0.7 21.75 0.64 0.84 0.6 0.7 19.7 0.6 1.2 

0.7 0.8 22.60 0.71 0.91 0.7 0.8 19.4 0.6 1.1 

0.8 0.9 22.28 0.74 0.96 0.8 0.9 19.8 0.6 1.1 

0.9 1.0 22.5 1.1 0.9 0.9 1.0 18.4 1.1 0.9 

1.0 1.2 17.1 1.5 1.5 1.0 1.2 17.00 0.40 0.94 

1.2 1.4 16.0 1.5 1.3 1.2 1.4 14.82 0.37 0.77 
Continued on next page 
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